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ABSTRACT 

Spectroscopic observations of distant quasars have resulted in the detection of molecular 
hydrogen in intervening damped Lyman-a absorption clouds (DLAs). We use observations 
compiled from different experimental groups to show that the molecular hydrogen abundance 
exhibits a dramatic increase over a cosmological time period corresponding to 13% to 24% 
of the age of the universe. We also tentatively show that the heavy element abundances in 
the same gas clouds exhibit a faster and more well-defined cosmological evolution compared 
to the general DLA population over the same time baseline. We argue that this latter point 
is unsurprising, because the general DLA population arises in a wide variety of galaxy types 
and environments, and thus a spans broad range of ISM gas-phases and abundances at the 
same cosmic time. DLAs exhibiting H2 absorption may therefore circumvent this problem, 
efficiently identifying a narrower class of objects, and provide a more sensitive probe of cos- 
mological chemical evolution. 

Key words: line: identification - ISM: molecules - Galaxy: abundances - intergalactic 
medium - quasars: absorption lines - ultraviolet: general 



1 INTRODUCTION 

The detection of cosmological chemical evolution is a cornerstone 
of modern cosmology, providing empirical details about galaxy for- 
mation and evolution, and independent supporting evidence of the 
big bang. High-precision measurements of chemical abundances at 
early cosmological epochs come from damped Lyman-Q absorp- 
tion systems (DLAs), arising in distant galaxies which intersect the 
lines of sight to high-redshift quasars. DLAs are easily detected in 
quasar spectra due to their particularly high neutral hydrogen col- 
umn densities [logj^Q A'^(H l) ^ 20.3 atoms cm^"^ by definition, 
though Voigt profile damping wings are seen at lower A'^(Hl)]. 
Generally, they comprise a clumpy distribution of clouds, spread 
typically over a few hundred kms^^. Detection levels for heavy 
element transitions in DLAs are ~7 orders of magnitude below 
that for A'^(Hl). High resolution optical spectra reveal transitions 
from many redshifted UV species, allowing detailed heavy ele- 
ment abundance patterns to be studied. The metallicity of DLAs 
- defined as the heavy element abundance with respect to hy- 
drogen, relative to that of the solar neighbourhood: [M/H] = 
logio[iV(M)/iV(H)] - logio[iV(M)/iV(H)]o - ranges over 2 or- 
ders of magnitude and is usually ~ 0.1 solar or lower. The presence 
of dust in DLAs is revealed by the relative depletion of refractory 
elements (e.g. Fe) with respect to those thought not to readily de- 
plete onto dust grains (e.g. Zn; Pettini et al. 1997), and by the sup- 
pression of blue light from the background quasar (e.g. Fall et al. 
Il989l) . These spectroscopic observations demonstrate that DLAs 
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arise along lines of sight through distant galaxies but they do not 
disclose details such as the galaxy's morphology, luminosity, mass 
or age. 

Evidence for an increase in DLA metallicity with cosmic time 
has emerged only g radually over the past few years f rom a number 
of different studies (Pettini et al. 1995; Lu et al. 'l99^:'Vladi loetalJ 
2OO0; Kulkarni & Fall 2002; Prochaska et al. 2003). The latter ref- 
erence provides the most compelling statistical evidence so far 
since, using the full sample of 125 DLAs, one sees a strong anti- 
correlation between [M/H] and absorption redshift, Zaba: Kendall's 
r is —0.32 with an associated probability of P(r) < 10~^. The 
need for such a large sample is due to the ~ 2 dex scatter in [M/H] 
at a given epoch. 

The types of galaxies responsible for DLAs can be inves- 
tigated using direct imaging methods or comparison of spectro- 
scopic data with models. The kinematics of optical resonance 
absorption transitions in low redshift DLAs appear consistent 
with 2-component models: a rotating disk plus an extended halo 
( lBriggsetalJll985h . A recent detailed radio study of one particu- 
lar DLA at Zshs = 0.437 clearly shows that at least some DLAs 
arise in the rotating disks of ordinary spiral galaxies iSriggs et alJ 
I2OOII) . High redshift (zabs ~ 3) DLAs have been interpreted simi- 
larly, using detailed kinematic modelling to argue that they are ro- 
tating ensembles of clouds with further evidence for this interpreta- 
tion coming from comparison of abundancej;citios in DLAs a nd the 
thick component of our Galaxy's disk jProchaska & Wolf9ll997L 
IWolfei200U . 

On the other hand, direct imaging of DLAs at Zabs < 1.5 re- 
veals the hosts to be a mix of irregulars, spirals and low surface- 
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brigh tness galaxies (LSBs; e.g. iLe Brun et alj 1 19971 : iRao et all 
l2003h . a result further borne out by a blind 21 -cm emission sur- 
vey at 2 = iRvan- Weber et al]l2003h . The observed number of 
DLAs per unit redshift interval and the A'^(Hl) distribution sug- 
gest that DLAs at Zabs < 2 are a roughly equal mix of LSBs and 
spirals; at higher redshift they are more likely to be dwarfs (that 
subsequently merge) since there are t oo few precursors of present- 
day d isk galaxies to explain the data I Boissier 2003: lBoissier et all 
l2003h . Fitting simple chemical evolution models to an observed 
A'^(H l)-[Zn/H] distribution suggests that at redshifts 2:abs ~ 2.5, 
about half of all DLAs are due to dwarf galaxies, while at Zabs < 1, 
giant galaxies may dominate the population ( Baker et al. 2000). H I 
21 -cm absorpt ion measurement s also hint that dwarfs could domi- 
nate at hie h-z iKanekar & Chengalur.20031) . 

Further work is clearly needed to generate an accurate inven- 
tory of DLA host galaxy types as a function of cosmic epoch. Nev- 
ertheless, the evidence for a broad mix of galaxy types is already 
compelling, if not conclusive. 



2 OUR CONJECTURE 

It is well established from studies of local galaxies that chemical 
abundances depend on morphology, luminosity and location within 
a galaxy. The diversity and evolution of the DLA host population 
will further complicate and distort estimates of chemical evolution 
over cosmological timescales. The DLA data themselves thus sug- 
gest that we may derive greater insight by attempting to disentangle 
these different galaxy types of DLA absorbers, rather than assum- 
ing all DLAs can be classified into one sample. A more sensitive 
test, or at least complementary test, for chemical evolution may be 
arrived at if we can select a subset of DLAs which arise in a narrow 
range of physical conditions. 

One possibility, which we suggest here, is to target those 
DLAs which exhibit molecular hydrogen absorption. There are cur- 
rently 8 examples of such H2 -bearing DLAs at redshifts Zabs > 
1.9. None of the absorbing galaxies have been identified directly 
through imaging, though hints at the galaxy types come from a 
comparison with our own Galaxy and its satellites. The mean 
molecular hydi'ogen fraction, / = 2Af(H2)/[2iV(H2) + iV(Hl)], 
seen along lines of sight through our the Galaxy is logjg / ~ — 1, 
and for the SMC and LMC it is log^Q / ~ -2. DLA detections 
lie in the range from — 6 < logj^o / < — 2 (upper limits on non- 
detec tions are —7 < logjQ / < — S. lLedoux. Petitiean & Srianand 
I2OO3I hereafter LPS03) but 2/3 of the det ections a re in the range 
-4 < login / < "2. As discussed in ILPSOBI > 90% of the 
sightlines through our Galaxy (and through the SMC) intercept 
H2 clouds. For the LMC, only 50% of sight-lines contain H2 at 
the same strength. DLAs exhibit H2 absorption ~ 20% of the time 
jLPS03h . and in this respect are more similar to the LMC. It is 
therefore possible that H2-bearing DLAs may arise predominantly 
in dwarf galaxies rather than massive spirals. However, given the 
rapid apparent molecular fraction evolution in Fig. 1, this interpre- 
tation is speculative. 



3 CURRENT OBSERVATIONAL EVIDENCE 
3.1 The data 

Table I summarizes the properties of the 8 H2 -bearing DLAs stud- 
ied here. H2 is detected in DLAs via the Lyman and Wemer-band 



UV lines and H I is detected via the Lyman series. Observed tran- 
sitions fall in the visible region for Zabs larger than ~ 1.9. Col- 
umn densities are measured from Voigt profile fits to absorption 
lines from each species. A compilation of results from H2 searches 
in DLAs is given in Table 8 of LPSOl for which 7 DLAs have 
confirmed H2 detections and metallicity measurements: 0013—004 
(z = 1.973), 0347-383, 0405-443 (2 = 2.595), 0528-250, 
0551-366, 1232+082 and 1444+014. The 2 = 1.9 73 DLA 
towards 0013—004 has been studied in detail bv .Petitiean et alj 
(2002): it comprises several absorbing components, of which 2 are 
dominant (their components c and d). These 2 components are dis- 
tinguishable in H2 but not in the saturated H I Lyman series. A'^(H l) 
has thus been measured at the mean redshift of these 2 components 
and we use the mean A'^(H2) in computing / for this DLA. The 
error w e use for / is derived from the range in A'^(H2) given in 
ILPSO3 !. We use the value and error for A'^(H2) in the DLA towards 
1232+082 from Srianand et al. (2000). 

We include one further DLA, 00 00-263. The H2 detection is 
regarded as only tentative in ILPS03I but has been carefully scru- 
tinised in lLevshakov et alj <200dl200ll) and relies on two H2 ab- 
sorption features, the L(4-0)R1 and W(2-0)Q(I) lines, the former 
appearing relatively free from Lyman-a forest blending. Since this 
is the highest redshift H2 -bearing DLA, the [M/H]-2abs and f-Zahs 
correlations derived below rely on this point relatively heavily. Fur- 
ther confirmation or refutation of this H2 detection is therefore de- 
sirable. A tentative H2 detection at 2 = 2.374 towards 0841 + 129 
has also been reported ( Petitiean et al. 2000), though confirmation 
requires future data and/or analyses. 

Here we omit the recent detection towards 0515— 4414 
iReimers et all2063) sinc e the local ionizing UV ba c kgrou nd may 
be relatively intense. See iMurphv. Curran & WebbI j^O^) for an 
analysis including this system. 

All metallicities are based on elements where depletion onto 
grains is expected to be small: 7 measurements are derived from 
[Zn/H], while the metallicity for 1232+082 is derived from [Si/H]. 
In most DLAs Si is observed to be slightly dust-depleted with 
respect to Zn so one might expect our conclusions below to be 
slightly strengthened if [Zn/H] could be used as a metallicity in- 
dicator for 1232+082. 

Fig. I illustrates the distribution of [M/H] and / for the 8 
DLAs with detected H2. The metallicity appears to change by 
~ 1.5 dex over the range 1.9 < 2abs < 3.4, while the molecular 
hydrogen fraction appears to change by ~ 6 dex. The scatter about 
the straight line fits is clearly larger than the statistical (1 u) er- 
rors in both panels. The validity of any simple parametric fit is thus 
dubious. Straight-line fits were nevertheless derived by quadrature 
addition of the same constant to each statistical error such that the 
normalised ~ ^- Least-squares fits give similar slopes and y- 
intercepts but our procedure results in reasonable error estimates 
for these quantities. 

3.2 Metallicity vs. Redshift 

With this procedure, a fit of the form [M/H] — As + B82abs to 
the 8 [M/H] points in Fig. 1 gives Ag = 1.0 ± 0.74 and Bg = 
—0.77 + 0.29. Using Kendall's r as a non-parametric correlation 
estimate gives r = -0.50 and P(r) = 0.08. That is, [M/H] and 
2abs are anti-correlated with 92% confidence. Com pare these statis- 
tics w ith those derived using the 60 DLAs from IProchaska et alJ 
(|2Qq3) in the redshift range 1.9 < 2abs < 3.4 which have Zn, O, 
S or Si metallicities: A = -0.75 + 0.41, B = -0.26 + 0.15, 
T = —0.13 and P(r) — 0.14. The slope, B, also compares well 
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Table 1. Summary of QSO and absorption line measurements for the 8 H2 -bearing DLAs considered h ere. Measurements of neutral hydrogen column density 
[N(Hl)], H2 column density [N(H2)], relative iron abundance [Fe/H], and metallicity [M/H] are from lLedoux et al] i2003h (LPS03) unless otherwise noted. 
The QSO magnitudes given are approximate V-band values. 'S/N (H2)' refers to approximate signal-to-noise ratio per ~ 2 km s~ ^ pixel for the continuum in 
the vicinity of the most clearly detected H2 absorption lines. 'M' refers to the element used to calculate the metalHcity. 
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17.9 


0347-383 


3.22 


17.3 
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2.595 20.90 + 0.10 

2.811 21.10 + 0.10 

1.962 20.50 + 0.08"= 

2.338 20.90 + 0.10^ 

2.087 20.07 + 0.07 
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+0.25 b 
0.06 
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0.10 
+0.23 
0.12 
+0.25 
0.19 



-5.71 
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+0.20 
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-3.82 

-1.48^ 



Notes: value for Q0528— 250 is not available in LPS03 but reasonably high S/N is expected from the QSO magnitude, exposure time and z^bs- 
References: ".Lu et al.. 1 1996>: 'Levshakov et al. ( 2000); ' Molaroet al. i,2OO0.) : ^ .Petitiean et al. I .200Z) : "=.Ledoux et al. ( 2002); ^ Srianand et aL i..200a) : 
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Figure 1. Upper panel: Metallicities for the 8 DLAs which exhibit H2 
absoiption (shaded circles), derived from the abundance of Zn (Si for 
1232+082). The relevant data are summarized in Table 1. The long-dashed 
line is a least-squares fit to the 8 points with expanded errors (see text). The 
dotted line is a similar fit excluding 0000—263. Small hollow diamonds il- 
lustrate m etallicity measuremen ts for the 60 DLAs between 1.9 < Zabs < 
3.4 from Prochaska et alj i2003i) . The short-dashed line is fit to the hollow 
diamonds, also with expanded errors. Lower panel: Logarithmic molecular 
hydrogen fraction. The dot-dash line represents the typical non-detection 
levels found bv Ledoux et al, l2Q03) . The age of the universe is calculated 
using = 0.30,nA = 0.70, Ho = 70kms~i Mpc"!. 



with that derived bv lProchaska et all for the overall 125 DLA sam- 
ple, i.e. -0.28 ± 0.05. 

Visually, the [M/H J-^abs distribu tion s uggest s two further re- 
sults: (i) As discussed in Ledoux et al. ( 2002, 2003), the H2-bearing 
DLAs have systematically higher [IVI/H] than the general DLA pop- 
ulation. A Kolmogorov-Smimov (KS) test reveals the two [IVI/H] 
distributions to have only 2% probability of being drawn from the 
same parent distribution. Removing the lowest and highest red- 
shift H2 -bearing DLAs, which are also the extreme [M/H] points. 
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Figure 2. As in the upper panel of Fig. 1 but for the distiibution of relative 
iron abundance. [Fe/H]. 



would give good agreement between the slope for the remaining 6 
points and the general DLA population, but the metallicity agree- 
ment would be somewhat worse (KS probability reduces to 0.8%). 
(ii) The []VI/H]-2abs correlation appears to be tighter for the H2- 
bearing DLAs compared to the general DLA population. To test 
this we subtract the linear fit to the data in each case and use an 
F-test to compare the variances of the residuals. The variance for 
the 60 general DLAs is 1.9 times that for the H2 -bearing DLAs, 
although the probability of this happening by chance is 39% given 
the sample sizes. 

In summary, the H2 -bearing DLAs display a faster increase in 
metallicity with cosmic time compared to the general DLA popu- 
lation. Furthermore, this increase is more well-defined in the H2- 
bearing sample and these DLAs seem to be chemically distinct 
from the general DLA population. 

These latter two results are also supported by the distribu- 
tion of relative iron abundances, [Fe/H], as shown in Fig. 2. Usin, 



[Fe/H] for the same 60 general DLAs from lProchaska et alJ 1200! 
as before, the KS probability for the two samples being drawn from 
the same parent distribution is only 0.01%. The F-test reveals that 
the variance of the general DLA population is 2.1 times that of 
the H2 -bearing sample but, again, this has reasonably low statisti- 
cal significance since we would expect a larger variance ratio 29% 
of the time by chance alone. Therefore, the H2-bearing DLAs do 
seem to form a chemically distinct sub-class of the DLA popula- 
tion which may encompass a narrower range of physical condi- 
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tions and environments. We note in passing tliat the slope of tlie 
[Fe/H]-2abs relation for the H2 -bearing DLAs, in contrast to the 
[M/H]-Zabs relation, is not significantly steeper than for the gen- 
eral DLA population. With the same fitting procedure above we 
find Ape = -0.78 ± 0.61 and the slope Bpc = -0.30 ± 0.24 
with Kendall's r = -0.21 and P(r) = 0.46 for the H2-bearing 
DLAs. The general DLA population gives A = —1.3 ± 0.36, 
B = -0.17±0.14, r = -0.10 and P(r) = 0.24. These shallower 
slopes are expected if dust depletion increases in importance with 
decreasing redshift, as implied by the /-Zabs relation discussed be- 
low (Section 3.3). We explicitly treat the evolution of the dust de- 
pletion factor in Murohv et al. (2004). 

Finally, it was earlier mentioned that the above conclusions 
rely relatively heavily on the highest redshift H2 detection to- 
wards 0000—263. If we remove this point, the straight-line fits to 
the [M/H]-Zabs relation give — 0.11 ± 0.81 and the slope 
Bj = —0.38 ± 0.33 with a non-parametric correlation estimate 
given by Kendall's r = —0.33 with P(r) = 0.29. Thus, the neg- 
ative [M/H]-z;abs correlation is not as steep, nor as robust when 
the 0000—263 data point is removed. However, note that the evi- 
dence supporting the hypothesis that the H2 -bearing DLAs form a 
chemically distinct population with a potentially narrower range of 
physical conditions is relatively independent of the 0000—263 H2 
detection. Comparing the 7 remaining H2 -bearin g DLAs with the 
49 gen eral DLAs between 1.9 < 2;abs < 3.1 frorn lProchaska et'al] 
t2Q03 ) gives a KS probability of 0.8% and the F-test gives a ratio 
of variances of 2.6 which occurs 23% of the time by chance alone. 

3.3 Molecular fraction vs. Redshift 

Turning to the lower panel of Fig. 1, we proceed as above and derive 
a straight-line fit for the /-^abs relation: Ai = 5.1 ± 2.4, Bf = 
-3.4 ± 0.9. Kendall's r is -0.50 with P(r) = 0.08. Thus, there is 
reasonable statistical evidence that the molecular hydrogen fraction 
increases dramatically with cosmic time. Previous studies do not 
explicitly identify this trend. In particular, although LPSOJ plot / 
versus z^hs for their sample of DLAs (their figure 16) they do not 
point out this apparently rapid evolution of / with Zabs- 

Again, removing 0000—263 from consideration weakens the 
statistical significance of the /-^abs correlation but the slope re- 
mains quite steep: the new fits with 0000—263 removed give Af = 
3.0 ± 3.0, B = -2.4 ± 1.2 with Kendall's r = -0.33 and 
P(r) = 0.29. However, we again stress the apparent reliability 
of the H2 de tection towards 0000—263 made bv lLevshakov et alj 
i200(l l200lh and point out that H2 systems with high-Zabs and 
high-/ (or, indeed, low-^abs and low-/) will need to be detected 
to show this correlation to be spurious. 

Below we consider possible interpretations of the f-z^hs trend 
but first discuss some systematic effects and biases. 



4 OBSERVATIONAL BIASES AND POSSIBLE 
SYSTEMATIC EFFECTS 

Could the strikingly steep apparent evolution of / with Zabs be sig- 
nificantly affected by observational selection effects? We consider 
some possibilities for both the [M/H]- and /-Zabs correlations be- 
low. 

Firstly, the H2 sample is inhomogeneous since the quasar 
spectra do not all have similar S/N and since the H2 detection meth- 
ods and criteria were not uniform. Indeed, the weak H2 lines de- 
tected towards 0000—263 and 0347—383 are at the same level as 



the many non-detections listed bv ILPSOsI (their figure 16). Typi- 
cally, these fall in the region illustrated in Fig. 1 by the light dot- 
dash line. The /-^abs correlation we report is therefore only ten- 
tative. However, it seems surprising that no high-Zabs/high-/ or 
low-Zabs/low-/ detections of H2 have been made. Secondly, the 
H2 detection limit will change with redshift: observed line widths 
increase as (1 + z), making higher-z lines more detectable, but 
the Lyman-a forest density increases roughly as (1 + z)^, making 
higher-2: lines harder to detect due to increased blending. This will 
be further complicated by higher order Lyman series forest lines for 
DLAs at low redshift compared to their corresponding QSO. Quan- 
tification of these competing effects requires numerical simulations 
which mimic the (unknown and, likely, non-uniform) techniques 
used to detect the 8 H2 systems studied in Fig. 1. 

The [M/H] and / values in Fig. 1 are not specific to the 
H2-bearing components because A*'(Hl) is determined from the 
damped Lyman-a line profile. The relative abundances of refrac- 
tory (e.g. Fe) and non-refractory (e.g. Si) elements are generally 
found to be uniform across the absorption profiles of most DLAs 
(Prochaska 2003), indicating that [M/H] is also likely to be uni- 
form. The notable exceptions to this rule are the H2 -bearing compo- 
nents where large dust depletion factors are found (e.g. LPS03 and 
the [Si/Fe] profile of 0347—383 in Prochaska 2003). These compo- 
nents usually dominate the non-refractory metal-line profiles and 
so, although the derived [M/H] and / measurements will be sys- 
tematically underestimated, the effect will not be large. Though the 
fitted slopes in Fig. 1 are likely to be reasonably robust against this 
effect, a larger sample and more detailed study of the metal-line 
and Lyman series profiles is required. 

Selection effects related to dust obscuration of quasars may 
also be present. In the H2 systems, / is somewhat correlated 
with dust depletion factor [Zn/Fe] >LPS03,) . providing evidence 
that the H2 is formed on dust grains, a v iew supported by the- 
ory jCazaux & T ielens"200l rGlovedl2003l cf. lLisztll2002h . How- 
ever, DLAs containing large amounts of dust (i.e. those with high 
/) could suppress detection of their backgroun d quasars a nd may 
therefore be 'missing' from our sample (.Fall et alJll98^ Since 
metallicity and dust depletion are also strongly correlated lLPS03h . 
this effect is likely to suppress, rather than create, the correlations in 
Fig. 1. Further to this, Boisse et al. ( 1998) and Prantzos & BoissieJ 
(2000) found an anti-correlation between [Zn/H] and A'^(Hl) in 
DLAs. For the 8 H2 systems studied here, there exists a mild cor- 
relation between A''(H l) and Zabs- If the latter is not simply fortu- 
itous, it should have resulted in a correlation between / and z^ha, 
opposite to that observed. 

The above implies that dust obscuration bias may somewhat 
suppress the observed /-Zabs correlation. However, several em- 
pirical studies suggest this bias is small. For example, DLA sur- 
veys using radio-selected quasars should not be subject to this 
bias. Recent work (E llison et al...2001.) indicates that the number 
of DLAs per unit redshift interval found in previous optical DLA 
surveys is not significantly biased, suggesting the number of 'miss- 
ing' DLAs is at most a factor of 2 below the radio-selected result. 
Nevertheless, we note 2 studies which tentatively suggest strong 
enhancements i n the incidence of abso rption towards radio-loud 
BL Lac objects jStocke&Rectol^ll997^ and optically red quasars 
tearilU et ajl998h . In the former case, many of the emission red- 
shifts are not known and so, if intrinsic absorption occurs more 
frequently in BL Lacs, the comparison with Ellison et al. (2001) is 
not useful. For the red quasars of iCarilli et aL 1,I99&') . only absorp- 
tion close to the qu asar redshifts was se arched for and so, again, 
the comparison with lElUson et alJ J200lh is tentative. Clearly these 
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possible discrepancies should be resolved observationally in future 
work. 



5 INTERPRETATIONS 

By selecting those DLAs exhibiting H2 absorption, one may fo- 
cus on systems with a narrower range of physical conditions than 
the general DLA population. Tentative support for this conjec- 
ture lies in the tighter metallicity-Zabs anti-correlation observed 
for the 8 H2 DLAs studied here. Hou et al. (2001) recently pre- 
sented detailed chemical evolution models which give a slope for 
the metallicity-Zabs relation of ~ —0.6 dex. They correct this pre- 
dicted slope for various observational biases to achieve the shal- 
lower slope observed for the general DLA population. However, 
the steep slope observed for 8 H2 systems may be less affected by 
these biases and may therefore be more comparable to the uncor- 
rected model slopes. H2-bearing DLAs could therefore provide an 
important probe of cosmological chemical evolution. 

That there exists such a large range (~ 6 dex) in the val- 
ues of / in Fig. 1 may not be surprising. For example, ISchavd 
i200 l') describes a photo-ionization model for clouds in local hy- 
drostatic equilibrium with a representative incident UV background 
flux, [M/H] and dust-to-metals ratio. The molecular fraction in this 
model shows a sudden increase of ~ 4 dex for only a small in- 
crease in the total hydrogen density. Therefore, the very steep /- 
Zabs correlation in Fig. 1 could be achieved with a modest increase 
in dust content at lower redshifts, consistent with the [M/H]-Zabs 
correlation. Within the Schaye model, one also expects an anti- 
correlation between the intensity of the UV background and /. 
However, the behaviour of the UV background flux with redshift 
over the range z — 2-3 is still a matter of considerable uncer- 
tainty. The strong decrease in / at high Zabs may also be con- 
sistent with recent Hi 21 -cm absorption measurements in DLAs 
teanekar & Chengalujiiool . where a generally higher spin tem- 
perature is found at z > 2.5. With an increased sample size and 
more detailed analyses, the f-z^bs correlation, if real, may provide 
complementary constraints on these problems. 



ACKNOWLEDGMENTS 

We enjoyed and benefitted from discussions with Ed Jenkins, 
Charley Lineweaver, Mark Whittle and Max Pettini. We thank 
the anonymous referee for critical comments which improved the 
manuscript. SJC gratefully acknowledges receipt of a UNSW NS 
Global Fellowship and MTM is grateful to PPARC for support at 
the loA under the observational rolling grant. This research has 
made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and 
Space Administration. 



Boissier S., 2003, in McWilliam A., Rauch M., eds, Carnegie Observato- 
ries Astrophysics Series Vol. 4. Origin and Evolution of the Elements. 
Carnegie Observatories, Pasadena, USA 
Boissier S., Peroux C, Pettini M., 2003, MNRAS, 338, 131 
Briggs F. H., de Bruyn A. G., Vermeulen R. C, 2001. A&A, 373, 1 13 
Briggs F. H., Wolfe A. M., Turn.shek D. A., Schaeffer J., 1985, ApJ, 293, 
387 

Carilli C. L., Menten K. M., Reid M. J., Rupen M. R, Yun M. S., 1998, ApJ, 
494, 175 

Cazaux S.. Tielens A. G. G. M., 2002, ApJ, 575, L29 

Ellison S. L., Yan L., Hook I. M., Pettini M., Wall J. V., Shaver P, 2001, 

A&A, 379, 393 
Fall S. M., Pei Y. C, McMahon R. G., 1989, ApJ, 341, L5 
Glover S. C. O., 2003, ApJ, 584, 331 
Hou J. L., Boissier S., Prantzos N., 2001, A&A, 370, 23 
Kanekar N., Chengalur J. N., 2003, A&A, 399, 857 
Kulkarni V. R, Fall S. M., 2002, ApJ, 580, 732 

Le Brun V., Bergeron J., Boisse P, Deharveng J. M., 1997, A&A, 321. 733 
Ledoux C, Petitjean R, Srianand R., 2003, MNRAS, 346, 209 
Ledoux C, Srianand R., Petitjean P, 2002, A&A, 392, 781 
Levshakov S. A., Kegel W. H., Agafonova 1. 1., 2001, A&A, 373, 836 
Levshakov S. A., Molaro P.. Centurion M., D'Odorico S., Bonifacio P., 

Vladilo G., 2000, A&A, 361, 803 
Levshakov S. A., Molaro P., Centurion M., D'Odorico S., Bonifacio P., 

Vladilo G., 2001, in Cristiani S., Renzini A., Williams R., eds, Deep 

Fields. Springer, Berlin, Germany, p. 334 
Liszt H., 2002, A&A, 389, 393 

Lu L., Sargent W. L. W., Bariow T. A., Churchill C. W., Vogt S. S., 1996, 
ApJS, 107, 475 

Molaro P., Bonifacio P., Centurion M.. D'Odorico S., Vladilo G.. Santin P., 

Di Marcantonio P, 2000, ApJ, 541, 54 
Murphy M. T., Curran S. J., Webb J. K., 2004, in Due P-A., Braine J., 

Brinks E., eds, lAU Symp. Ser Vol. 217, Recycling intergalactic and 

interstellar matter. Astron. Soc. Pac, San Francisco, CA, U.S. A 
Petitjean P, Srianand R., Ledoux C, 2000, A&A, 364, L26 
Petitjean P, Srianand R., Ledoux C, 2002, MNRAS, 332, 383 
Pettini M., King D. L., Smith L. J., Hunstead R. W., 1995, in Meylan G., 

ed., QSO Absorption Lines. Springer- Verlag, Berlin, Germany, p. 71 
Pettini M., King D. L., Smith L. J., Hunstead R. W., 1997, ApJ, 478, 536 
Prantzos N., Boissier S., 2000, MNRAS, 315, 82 
Prochaska J. X., 2003, ApJ, 582, 49 

Prochaska J. X., Gawiser E., Wolfe A. M., Castro S., Djorgovski S. G., 

2003, ApJ, 595, L9 
Prochaska J. X., Wolfe A. M., 1997, ApJ, 487, 73 

Rao S. M., Nestor D. B., Turnshek D. A., Lane W. M., Monier E. M., Berg- 
eron J., 2003. ApJ, 595, 94 
Reimers D., Baade R., Quast R.. Lev.shakov S. A., 2003, A&A, 410, 785 
Ryan- Weber E. V, Webster R. L., Staveley-Smith L., 2003, MNRAS, 343. 
1195 

Schaye J., 2001, ApJ, 562, L95 

Srianand R., Petitjean R, Ledoux C, 2000, Nat. 408, 931 
Stocke J. T., Rector T. A., 1997, ApJ, 489, L17 

Vladilo G., Bonifacio R, Centurion M., Molaro R, 2000, ApJ, 543, 24 
Wolfe A. M., 2001, in Funes Jose G. S. J., Corsini E. M., eds, ASP Confer- 
ence Series Vol. 230, Galaxy Disks and Disk Galaxies. Astron. Soc. Pac, 
San Francisco, CA, U.S.A, p. 619 

This paper has been typeset from a TgX/ BTgX file prepared by the author. 



REFERENCES 

Baker A. C, MathUn G. P, Churches D. K., Edmunds M. G., 2000, in Fa- 
vata F, Kaas A., Wilson A., eds, ESA SP Vol. 445, Star Formation from 
the Small to the Large Scale. European Space Agency, Noordwijk, The 
Netherlands, p. 21 

Boisse P, Le Brun V., Bergeron J., Deharveng J., 1998, A&A, 333, 841 



